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Consider a complex energy z for a Ai'-particle Hamiltonian H and let x be any wave packet 
accounting for any channel flux. The time independent mean field (TIMF) approximation of the 
inhomogeneous, linear equation {z — -ff)|'I') = |x) consists in replacing ^' by a product or Slater 
' determinant (j> of single particle states ipi. This results, under the Schwinger variational principle, 

, into self consistent TIMF equations [r]i — hi)\tfi) = \xi) in single particle space. The method is a 

, ^ ■ generalization of the Hartree-Fock (HF) replacement of the A^'-body homogeneous linear equation 

^ (E — H)\'^) = by single particle HF diagonalizations (e^ — hi)\ipi) — 0. We show how, despite 



00 



Oh' 



strong nonlinearities in this mean field method, threshold singularities of the inhomogeneous TIMF 
equations are linked to solutions of the homogeneous HF equations. 



I. INTRODUCTION 



fT^ I After the success of the mean field approach for bound state systems in various fields of physics, it was only 
natural to try the mean field concept for scattering states as well. The original attempt [1] was the time-dependent 
Hartree-Fock (TDHF) method where one solves the single-particle equations of motion as initial value problem in 
time. From the resulting solutions at various impact parameters one may then calculate the classical cross section. 
With no specification of the final state, the method is restricted to inclusive reactions. A serious, conceptual problem 
arises from spurious cross channel correlations [2,3]: when projecting the TDHF Slater determinant for large times on 
an orthogonal set of channel wave functions, the expansion coefficients and the respective S- matrix vary in time ad 
infinitum. To overcome the shortcomings of TDHF, the time-dependent mean field (TDMF) approach [2,3,4] expands 
' the density in two sets of (bi-orthogonal) single-particle wave functions and solves the equations of motion as boundary 
• value problem in time, fixing initial and final densities. It has been proven that for TDMF an S-matrix can be defined 
' which becomes asymptotically constant [2]. The problem with TDMF lies in combining self-consistency with given 
• • • boundary conditions in time [3,5]. No practicable algorithm for this highly "non-local" problem exists upto date for 
. !^ I use in actual numerical calculations. A third approach is the time-independent mean field (TIMF) method [6], based 
' on a Schwinger-type variational principle [7] for matrix elements of the resolvent or T-operator between given initial 
H and final states. The method uses two sets of variational single- particle functions, analogous to TDMF, and leads to 
- - 1 inhomogeneous equations of Hartree-Fock type which can be solved iteratively for given total energy of the system. 
TIMF is free of the conceptual and practical problems of TDHF and TDMF, resp., and has been tested successfully 
on a number of simple systems. It can be extended to incorporate particle-hole correlations, as has also been done for 
TDHF, within a generalized random-phase-approximation [8]. The present paper goes beyond the above problems [9] 
by studying the continuum singularities of this TIMF approach for collisions. 

Consider a finite number N of particles. Factorized wave packets (shifted Gaussians in momentum representation 
for example) make an overcomplete basis in their Hilbert space of wavefunctions. Hence the calculation of a retarded 
Green's function amplitude T) = {x\{z — H)^^\x), where z) x is a product, \x) — YiiLi IXi)? s-nd each single particle 
wave function Xi is real rather than complex, makes a fully generic problem. Such factorization simplifications are 
not physically restrictive and help in the analysis of a mean field theory of collisions, the subject of this paper. 

For the sake of simplicity, we deal yet with spinless, distinct particles only and short range interactions Vij for the 
Hamiltonian 77 = X]iLiPi/(2"^i)+Si>j=i ''^y ■ The case ofidentical particles can be treated later and, in the following, 
any reference to a Hartree method may be understood as a reference to a Hartree-Fock (HF) method if necessary. 
Again for simplicity, we consider the calculation of diagonal collision amplitudes only, {k\V{E^ — H)^^V\k), Born 
term subtracted. Generalizations to distinct prior and post interactions, V,V', are kept for future work. The state 
\k) is taken as a plane wave of relative motion in any two cluster channel ground state and the product \x) = V\k) is 
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a square intcgrablc state in the A^-particlc space. Finally z is any complex number E + iF, and the usual limit S+ at 
the end of any calculation reads F +0. 

It is trivial to use the Schwinger variational principle [7] and show that V is the stationary value of the functional 

under variations of The corresponding Euler-Lagrange equations read, with retarded boundary conditions and 

arbitrary norms and phases of and , 

{z-nm = \x). (*'|(^-if) = (x|. (2) 

The variational equations which occur in the time independent mean field (TIMF) [6] theory of collisions read, 

{m-K)Wi) = \xi), - K) = {xrV i = i,-,N. (3) 

They are obtained from Eq. (1) when x, and the approximation cj), resp. (j)' , chosen for 'f, resp. are products of 
single particle orbitals, XiTfi^fiJ respectively. Such TIMF equations are very simple [6]. Except for a single particle 
density operator p defined non diagonally as p{f',f) = '^j^'Pi{r')(p'*{r), they arc just Hartree(-Fock) equations 
completed by a right hand side, representing the image of the channel in single particle space. In the following, is 
restricted to such factorized source functions x and trial functions (!>, (f)' and will be labeled F. A saddle value under 
such a restriction of 0' is not necessarily unique any more. It will be denoted by D instead of T) and may request an 
additional, identifying label. Now our claim is: bound and unbound solutions of the usual Hartree(-Fock) equations, 

{ei - hi)\<pi) = 0, (4) 

induce singularities of the one-body variational conditions, Eqs. (3). 

This reminds, naturally, of the strict connection between the singularities of the linear, inhomogeneous problem 
{z — H)\^) = \x) in the A''-body space and the solutions of the linear, homogeneous Schrodinger equation {E—H)\^) = 
in the same space. Because of the nonlinear nature of Eqs. (3-4) in single particle space, our claim is not obvious, 
and will be qualified in this paper. 

Actually, in a previous paper [10], the claim was already substantiated in part : those energies Eh, for which a 
bound Hartree(-Fock) solution (t)H is found, generate poles of the approximate amplitude D provided by saddle points 
of the restriction F. Furthermore {z — Eh)D \{x\(I^h)\^ / {(I>h\<I>h) when z Eh- Despite the nonlinearity of the 
approximation, such a residue at such a pole is almost expected. The analogy with the poles of V at exact eigenvalues 
for bound states is striking. We are now interested in a more difficult question, namely, is there a similar analogy at 
higher energies, when singularities of scattering and rearrangement collisions (thresholds, cuts) occur? 

In Section II we briefly recall a very simple, soluble model [11], used earlier among several other models to validate 
D as an approximation of V. The model is reintroduced for pedagogical reasons first, to illustrate a derivation of 
Eqs. (3). Then, and mainly, it is used to provide a complete investigation of singularities, for it boils down to 
manipulations of polynomials. In Section III we introduce an enriched model, exactly soluble too. Section IV contains 
a generalization and discussion of the results obtained in Sections II and III. Finally Section V contains our conclusion. 



II. FIRST MODEL, BARE PROPAGATION, SYMMETRIC MEAN FIELD, TWO-BODY THRESHOLD 

In this soluble model, there are only two one-dimensional particles with just their kinetic energies, and different 
masses rUi = l/(2aj), hence H = aip\ + 02^2- While the inversion of 2 — is numerically trivial and allows a good 
validation [11] of the TIMF approximation, the formal expression of {z — H)~^ in terms of one-body propagators 
{r]i — aip\)~^ and (772 — ^2^2)"^ is less trivial, as it demands a convolution. The TIMF method consists in replacing 
the convolution by just one product, namely 

{z - aipl- a2pl)~'^\XiX2) oc {m - a.iPi)~^\Xi) - a.2pl)~^\X2) ■ (5) 

This comes from variations 5/6(pi of the functional F. An additional simplification results from a further remark : in 
those representations where x and H are real, one finds from Eqs. (2) that j^*') = j^**), hence the possibility of just 
one trial function ^ if one uses a Euclidian ( | ) rather than a Hermitian ( | ) metric, 
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- (vi/|(z-F)i*) mz-Hm' 

For the present two particle model, the factorization of x into two single particle wave packets with real wave functions 
Xi , X2 allows us to use the following form of F, 



{(P1(P2\ {Z - aipl - a2P^) \Vl'P2) 2:((^l|(^l)(<^2|</32) - (<^l|aipf |<pl)((/J2|</'2) - {ip\\ipi){ip2\aiP2\V2) ' 

We assume that Xi)X2 are real in the momentum representation. The functional being insensitive to the norms and 
global phases of , <^2 , elementary manipulations of 5F/ 6ipi yield, in the same momentum representation, 

'Pi{p)= i=l,2, (8) 

Vi - 0.iP 

with 

J dpyj{p)ajp^ J dpx^[p)[ajj>^ -ijj)-^ J ■ 10 A -^^ 

Vi = z jr- — — = z-rjj- -TT-j — 27T7 2 W2 ' * = 1' 2, and j = 1, 2, and j ^ i. (9) 

It is convenient at this stage to define the integrals, 

Ii = -{XM)= [ dp ^l^^ , i = l,2, (10) 
J ttip^ - rii 



and notice that Eqs. (9) then read. 



^i^=--E^^' i = 1.2. (11) 

^ i=l 

If furthermore one defines auxiliary variables uoi by the conditions 

r?i = ajW?, Swj > 0, « = 1, 2, (12) 
then it is useful to define Jj = ajlj. And Eqs. (11) become, 

'^aju;jJj-^ = z-J2<^i'^l i = l,2, (13) 
where a contour in the upper half plane of the complex variable p defines the integrals, 

The special cases Qijj while ^rij > define cuts in the complex Tjj-plane. These correspond to Qcoj ^ in the 
Wj- -plane. 

When the two particles are identical, it may be interesting to symmetrize and antisymmetrize Eqs. (13) as, 

2 2 
'^aiU}iJi—^ = z-'^aiU)f, (15) 
i=i ^-^^ i=i 

and 

ai wi Ji — - - a2a;2 J2 -Ti- = 0, (16) 
dJi (IJ2 

and identify cases where the mean field might break their symmetry. But we shall keep the particles, and/or their 
channel wave packets, distinct for a while. 



3 



A soluble model, involving only the manipulation of polynomials, is obtained if one chooses the forms of the wave 
packets as follows. 



Xiip) 



7j 



1/2 



1,2, 



yielding the simple result. 



(17) 



(18) 



Resulting polynomial equations turn out to have a lower degree if Kj = 0, for then Jj becomes Jj = —l/[LLij{ijjj 
As will be found in this and the next Sections, two kinds of singularities emerge, t) "physical" ones, which essentially 
depend on z and are not very sensitive to "technical" parameters Kj,aj,'yj, and ii) "technical" singularities, more 
sensitive to such parameters. The analytical continuation provided across 77 cuts [12] by this uj representation is clear. 
Once Eqs. (13) have been solved, the saddle point values of the functional read, using Eqs. (7-11), 



D 



{aiLof + a2w| - -z) Ji J2 



aia2 



(19) 



The search for singularities of -D as a function of the physical energy z thus consists in eliminating oji , LU2 between 
Eqs. (13) and Eq. (19). The former read, after elementary manipulations which take advantage of Eq. (18) when 
Ki=K2 = 0, 



2aix'^y + aij2X^ - 0272?/^ + 2y2 + 72-2 = 0, 2a2y'^x + 02711/^ - ai7ia;^ + 2xz + jiz 



0, 



(20) 



where it was convenient to set u)i = ix, ^x > 0, and C02 = iy, ^y > 0. Equivalently, if we scale x and y into x = 71a;' 
and y = 72?/', respectively, the same equations read, 



{Aix'^ + z){l + 2y') - A2y'^ = 0, (Ast/'' + ^) (1 + 2a;') - Aix'^ = , 
with Ai = oi7i ^iid A2 = 0272- elimination of y between Eqs. (20) gives, 
2af7ia;'' — af{4:Z — ai^f + 0272)2;^ — aiz{8z — ai^f + 40272)2;* — 

2ai7i^;(32: + a27|)a;^ - z'^{4:Z + ai'yf + 4a27|)a;^ - 4^;^7i(^; + a27|)a; - 'yiz'^{z + 027!) = 0, 

72 (ai + 2x z + ji z) 



y= - 



(2.x + 71) (ai.i;2 + z) 



(21) 



(22a) 
(22b) 



The multiplicity of solutions is thus 7, which raises a problem for the identification of a solution, if possible unique, 
which accounts for a physical approximation. 

In turn, if one inserts Eq. (18) into Eq. (19), one obtains D as a rational function of a;i,ix;2, or x,y as well. Upon 
taking advantage of Eq. (22b), this rational fraction reduces into a rational fraction of x only, hence a polynomial 
relation between D and x, with degree 1 for D, 

a^a272a;^(7i + xf{aix^ + 71Z + 2xz)D = Aalx'^{x + 71) + afiai^^ + 0273 + ^'2z)x^ + Uaf^izx'' + aiz{3aijf + 
4 027I + 122:)x* + 2oi7iz(a272 + 6z)x^ + 2^(3ai7i + 4a27| + 42)a;^ + 471^^(0272 + z)x + 71^^(0272 + z). (23) 

The same result is obtained if one uses Eq. (7) instead of Eq. (19). 

An elimination of x between Eq. (23) and Eq. (22a) finally gives a direct, polynomial condition relating D and z, 

zlzl{zi + 1)(^;2 + 1){Z1 + Z2 + - ^ZiZ2{zx + 1)(2;2 + l)(-2l-22 - 4^1 - 4^;2 - 4)56 _ 4[32;3^2 ^ 3_j,2_2,3 ^ 

(20^^2:2 + ft^zlzl + 2Qzizl) + {l&zl + S,S,zlz2 + 88^1 2;| + 16^) + (32^^ + 84^i2;2 + 32^1) + 16(^1 + Z2)]D^ + 
16[3zfz^ + 39(z?Z2 + ziz|) + (32z? + 91zi22 + 32z|) + 48(zi + za) + 16]D^ + 



16[3ziZ2 + iziA 



izf + 91ziZ2 + 8z|) - 88(zi + Z2) - 64]L)3 + 



192[-ziZ2 + 4(zi + Z2) + 8]D^ ~ 64(zi + Z2 + 16)D + 256 = 0, 



(24) 



where Zi = ai-yf/z and D is scaled as D = D/z. The degree 7 for x in Eq. (22a) is correctly refiected here by the 
same degree for D (and D). Conversely, given an amplitude D, the degree of the polynomial condition, Eq. (24), 
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with respect to z is 4. Hence there are 7 approximate amplitudes offered by TIMF for each energy, while the inverse 
problem, "given the TIMF amplitude, find the energy" , has 4 solutions. 

This model, although soluble, thus creates a complicated Riemann surface. Criteria are necessary to select one 
physical sheet, or physical pieces of sheets. Obvious candidates are the conditions > 0, Sfty > when Eqs. (20) are 
solved. Concerning Eq. (24), the very definition of V demands that T) be real and negative if z is real and negative. 
When 2: > with a slight and positive imaginary part, then SjD must be negative. Those roots D which show the 
same properties should thus help the identification of suitable sheets. 

The argument is made much simpler if the "technical" parameters a\^\ and 027! taken equal to some common 
value Q. This amounts, in some sense, to consider identical particles, although a\ may still differ from 02- Then Eq. 
(24) factorizes as 

[4 - 4(61 + z)D + e{e + z)D'^\ ^ [16 + 8(36» - ^z)D + 4(361^ + IWz + ^z^)D'^ + 6i^(26» + z)D^\ = 0. (25) 

If is used as a unit for z and similarly 1/6* is used as a unit for D this reads as well, 

[4 - 4(1 + z)D + (1 + z)!?^] =0, (26a) 
[16 + 8(3 - 4z)D + 4(3 + lOz + '^z^)D'^ + (2 + z)D^\ = 0. (26b) 

The presence of a squared polynomial as the first factor in Eq. (25) reflects a "symmetry breaking" by the mean field 
approximation. Indeed, when analyzing the corresponding solutions of Eqs. (20), one finds that each pair of roots 
with X ^ y. is accompanied by a pair generating the same value of D. Such a degeneracy thus makes, 

out of 4 of all the 7 solutions for {x,y}, two distinct values for D. All told, D then takes 5 distinct values. The 
remaining 3 solutions account for the degree 3 present in the second factor of Eq. (25). It is easy to verify that such 
3 solutions are "symmetric" , namely x = y. Notice that the symmetry breaking generates a rational inverse function, 

_ (Afc-2)^ .27^ 

while the symmetry conservation generates an equation of degree 2 for z. Since Zbk must be counted twice, one recovers 
the 4 solutions of the inverse problem. 

It turns out that the "symmetry breaking" sector violates the double condition, 5R.t > 0, SRy > 0. Hence the 
properties of this sector are listed in an Appendix only. Turning now to the symmetric amplitude Dgy, the choice of 
a physical branch is reasonably easy, see Figs. 1-2. In Fig. 1, the lower half plane contains a loop acceptable as a 
physical candidate. We verified that > for this loop. Despite a suitable ^x > if 2 > 0, the other branch in 
Fig. 1 is clearly not acceptable, for it contains values Dgy with positive imaginary parts. Nor can one accept the third 
branch, seen in Fig. 2, despite its correct sign for ^Dgy For it violates both the limit I? when ^ 00 and the 
obvious condition "3?I? < if .2 is real and negative" . Furthermore SRa; is found unsatisfactory for this third branch. 
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FIG. 1. Complex D-plane. Two trajectories of symmetry conserving amplitudes as functions of SRz when — 1. Growing 
blue dots: 5Rz grows from — oo to —0. Growing red dots: increases from +0. The third trajectory lies far in the lower half. 
Color is available online at www-spht.cea.fr/articles/t02/148/ 

Two values of z generate branching for Dgy. With a single root D = —1/5, reasonable, a double root D = 16, 
unphysical, occurs for z = -27/16, with expansion D = 16 - 256/9 (z + 27/16) ± 8192/243 [-{z + 27/16)3]^/^ Hence 
a familiar square root cut can be used to disentangle the two corresponding sheets, both unphysical. The value 
z = —27/16 docs not represent a natural threshold for the present model. More physical, obviously, is the triple root 
singularity, D = —2, which occurs at z = 0, the true threshold. It is illustrated by Fig. 3, where a tiny imaginary 
part 3z = .0001 was added in order to separate branches. It will be noticed here that, although the physical branch 
gives real values of Dsy when z < and complex values of the same when z > 0, there are always one real root and 
two complex conjugate roots on both sides in the vicinity of z = 0. This happens indeed because the corresponding 
discriminant, A3 = 256z^(27 + 16z)'^/(2 + z)*, actually changes sign, not for z = 0, but rather for z = —27/16. This 
helps to understand the nature of the unphysical singularity occurring at z = —27/16. It gives an early "warning" 
of the (cubic) physical threshold singularity, z = 0. An elementary, but slightly tedious calculation provides the 
expansions of the 3 branches in the vicinity of z = 0, namely D = —2 — j 2^/"^ 3 z^^^ + ©(z^/"^), where j is either 1, or 
any one of its complex cubic roots (—1 ± i Vs) /2. 
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FIG. 2. Complex D-plane. All trajectories Day when '^z = .075. Scales of trajectories made compatible by replacing radii 
from the origin by their square roots. Hence, for instance, announcing the double root D = 16 when z = —27/16, blue branches 
cross each other near yfD — 4. Growing blue dots: grows from — oo to —0. Growing red dots; SRz increases from +0. 
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FIG. 3. Complex D-plane. Triple merging Day —2, with Qz = .0001, —.09 < SRz < .09. Lower right branch physical. 



7 



Consider Eqs. (21) and set Ai = to factorize the resultant, Eq. (22a). Then scale y and z as proportional 
to 7i, 72 and Ai, respectively. For the sake of simple numbers, this strictly amounts to set a common value a\ = 
a2 = 7i = 72 = 1, hence A\ = A2 = 1, for those reduced equations which govern the scaled variables and parameters. 
For the symmetry sector, x = y, both equations, Eqs. (20), then boil down to 2x^ + 2xz + z = 0. It is trivial to 
find that, at threshold z ^ 0, all three roots have a leading term x = + 0{z'^/^), while, as already found, 

Dgy = — 2 + 0{z^/^). Obviously, below threshold, one must select the real root x, which gives a real amplitude. 
Conversely, above threshold, one must select that complex x which gives a retarded amphtudc. 

All told, for Dgy, cuts needed in the z-plane are a cut from to +00 for the cubic branching and, for instance, a 
"technical" cut from —00 to —27/16 to create an additional seam between the second and the third sheets. 

Now we consider additional cuts, namely those created by the condition 3cj = 0, or, identically, by the condition 
3?a; = 0. These occur because the solutions of realistic problems demand numerical, iterative calculations of rii and 
ipi before obtaining D. This means inversions of operators (ajW? — hi) in sequences of successive approximations of 
w's (and self consistent /I's when potentials arc involved). Obviously, every time vanishes or becomes too small, 
numerical precautions are in order. Also, since the physical energy is on shell, z = E + iO"'", with a retardation 
boundary condition for many-body propagation, one would feel more comfortable with retardation also for the single 
particle energies rj oc u'^. Advanced ij's are not to be ruled out a priori, because it is well known that mean field 
approximations can be excellent while breaking many-body symmetries. But, clearly, branches of x's which cross such 
cuts 5Ra; = deserve some cautious scrutiny. 

For the present case where Ai = A2 for "symmetric" bare propagations, and still with simple numbers ai = ji = 1, 
our results are shown in Figs. 4-5. (For the academic, "symmetry breaking" case, see the Appendix with Figs. 13-14.) 
Fig. 4 is a contour plot of the product ^xi^X2^X3 of the real parts of the 3 roots as functions of z in the z-plane. 
Darker areas indicate an increasing positive product (two out of the three 3?x's are < 0), while the lighter areas mean 
a more and more negative one (one negative 3?a; only). The product vanishes along the contour line separating the 
light grey area from the moderate grey one. It will be noticed that this line contains the point z = 0. Hence the cut 
relevant to D and that relevant to x's both contain the two-body threshold. Notice, however, that, except at such a 
treshold, a real z induces complex r^'s. Namely, propagation energy cuts do not follow the real axis in the z-plane. 
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FIG. 4. Complex z-plane. Cut caused by the condition Ka; = for symmetry conserving roots. The cut is the contour line 
separating the lighter grey area from the darker grey one. 

The next Figure, Fig. 5, shows the trajectories of the roots when we freeze 5iz = 0.1, above threshold, and let ^z 
run from —1 to +1, hence allowing one ^x, then a second one, to change their signs. The sizes of dots are coded as 
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follows: minimal for SJz = — 1, growing until 3z = 0, minimal again for small positive values of Qz, then growing again 
until 3z = 1. The lower branch is the best candidate for physical roots, because it provides a growing retardation, 
< 3?7 = 3(— x^), when is positive and grows. As predicted from Fig. 4, there is an interval for where 2 roots 
X have a positive ^x. 




FIG. 5. Complex a;-plane. Trajectories of the symmetry conserving 3;'s when SRz — 0.1, while Qz crosses the cut shown by 
Fig. 4. Blue dots growing when grows from —1 to 0. Red dots growing when Qz grows from to 1. 

To conclude this Section, the main result derived from this elementary model with bare propagation of two particles 
lies in the systematic, physical, two-body threshold found at z = in the energy plane (z-plane) for all the mean field 
quantities, whether amplitudes D or propagation energies 77. This threshold is, obviously, a common feature of both 
the exact problem and the corresponding Hartree problem. For amplitudes D, a cut in the z-plane extends from the 
threshold to +00, as seen in both the "symmetric" and "breaking" submodels. For propagation energies 77, the cut 
starts from 2 = 0, indeed, but deviates from the real semi-axis. For both D's and rj's, the cost of the nonlinearity of 
the TIMF approach is reflected in additional, unphysical, "technical" singularities. But such unphysical singularities 
are not beyond interpretation either, as shown by the analytical properties listed in this Section. Incidentally, as 
discussed earlier [13], unphysical singularities may be washed out by a linear admixture of the various solutions of the 
nonlinear mean field problem. The next Sections will show even better how physical cuts remain a significant feature 
of the TIMF approximation. 



III. SECOND SOLUBLE MODEL, ONE-BODY THRESHOLD 

Here again we consider two one-dimensional particles, and particle 2 is still free with a pure kinetic energy /12 — CL2P2 
for its Hamiltonian. But now the complete Hamiltonian H = hi + h2, while still separable, involves a bound state for 
particle 1, because we set hi = aipf — A|xi)(xil: with an attractive enough potential. For technical reasons which will 
soon become clear, the form A|xi)(xi| of this potential makes use of the same wave packet Xi taken as a channel wave 
packet. The numerical inversion of z — iJ is still easy and allows another good validation of the TIMF approximation. 
The formal expression of (z — H)~^ in terms of one-body propagators {r]i — hi)~^ and (772 — ^2)^^ demands again a 
convolution and the TIMF method consists in replacing the convolution by a product, 

{z - hi - h2)-'\xiX2) (X ivi-hi)-'\xi)im-h2)-'\x2). (28) 
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This comes again from variations S/Sipi of the functional F. And a further remark can be repeated: in those repre- 
sentations where x and H real, we obtain = see Eqs. (2). Hence the possibility of just one trial function 
(f) under a Euclidian rather than a Hermitian metric, see Eq. (6). The factorization of x into two real wave packets 
Xi,X2 essentially retains Eq. (7), which actually becomes, 

p ^ {XlX2WW2f ^ (Xi|y^i)^(X2|y^2)^ .29) 

{^pl<f2\ {z-hi- /l2) Z{ipi\ipi){ip2\ip2) - {^l\hi\ipi){ip2\^p2) - {Vl\V>l){V>2\h2\ip2)' 

We use the same xi,X2, real in the momentum representation. The fimctional being always insensitive to the norms 
and global phases of tpi, (^2, the same manipulations of dF/S(pi yield, in the same momentum representation, 

IVi) = im - /ii)-'lXi), iVi - aipDl^i) = Ixi) - A|xi)(Xi|Vi), Mp) = [1 - A(xibi)] , (30) 

r]i — ttip 

'P2[P) = (31) 

where it is better, temporarily at least, to retain the factor 1/ = [1 — A(xi|vi)] for ifi. The same quantity v, as will be 
seen shortly, cannot be discarded from the self consistency conditions of the pair r]i,r]2, 

JdpipUp) a2P^ = ^ _ ^2 - ■f'^P^2ip)ia2p'^-V2)~'^ ,32) 
^ Jdpiplip) ^ J d.pxl{p){a2p'^ -V2)~'^' 

{ipi\hi\ipi) ((^i|(/ti - r]i)\(pi) (Xil(fei -Vi)~^\Xi) ioo\ 
V2=z — ■ — —=z-r]-i — ■ — =z-'ni-- — — — — -. {66} 

(Vim) mm) (xiK/n -m) Ixi) 

Indeed, it is necessary to consider the matrix clement, 

2:l = (Xl|(/^l-^l)"'IXl), (34) 

and notice that Eqs. (32-33) become, 

r dv2 iX2\^2) dr]i (Xilvi) .oK^ 

dh (1^2 1 1^2) all mm) 

The integrals /i,/2 were already defined by Eq. (10). Returning to Ji, and to the factor v which accounts for the 
separable potential present in h\, an elementary manipulation of Eq. (30) gives, 

(Xibi) = -Xi, I, = ^Jl_. (36) 

Again we define auxiliary variables iVi by Eq. (12) and integrals Jj = Ujlj by contours in the upper half plane of the 
complex variable p. Then Eqs. (35) become 

d(jJ2 2 2 r. ^ CLjl o 2 

2 02^2^2 -7^ = z — aicoi — a2U)2, 2aiWiX\ —— = z — aiw^ — 02^2. (37) 
dJ2 dl\ 

It will be recalled here that a (unique) bound state occurs for hi for any positive value of A, at an energy % < 0) 
defined by the well known condition. 



\=ldp^^='-l^^ 
A J aip'' - rjo ai 



dp 2 = ' % = "li^o' ^'^o = 0, Swo > 0. (38) 



Indeed, the right hand side is monotonically increasing when rjQ nms from —00 to and the same r.h.s. diverges at 
rjo = 0, see Eq. (18), because of our choice of a Lorentzian form for Xi- Accordingly, an explicit form of Eq. (38) is, 

A + aia;o(wo + «7i) = 0, SRwo = 0, ^u>o > 0, (39) 

or, in terms of r]o, 

(r/o + A)2 + 0172 r]o = 0, (40) 
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with obvious scaling properties. (Indeed, if the scale is set by A for instance, it is convenient to define Ai = aijf, 
and the relevant scales are, obviously, Ai/X and r/o/A.) Threshold singularities are expected for Eqs. (32-33) when z 
reaches the one-body threshold r]o, besides the already found two-body threshold z = 0. 

The saddle point value D deduced from Eq. (29) reads, upon taking advantage of Eqs. (30-37), 

D = {aiul + a-iJi - z) Ji h ■ (41) 

This formula, Eq. (41), is an obvious generalisation of Eq. (19). In the same way as we did in the previous Section, 

we shall again eliminate ui and lu2, or rather the strictly equivalent variables x = —icui and y = —iuj2, between Eqs. 
(37) and Eq. (41). It is then useful to define a parameter A2 — 027!, quite similar to Ai and it is also easy to predict 
that the solution D{z) scales in terms of Ai/X, A2/X, z/X and XD. The Lorentzian choice for Xi,X2, induces the 
following forms for Eqs. (37), when we replace 01,02 by ^1/71,^42/72, respectively, 

Ai7|x2 + 2Ai-f2x'^y - ^27??/ + ihh + '^lh2yz = 0, 
2717IAX + A27iy2 + 2A27ixy2 - Ai'^jx^ + -il^i^z + 27i7|a;z = 0, 
y ^ {X-il+Aix'^)x + {-ii+2x)-ilz 
72 (71 -F2x)(^ix2 +7^2) 

These scale obviously in terms of a;/7i and y/^2- It is then convenient to set 71 = 72 = 1 in Eqs. (42). 

Simultaneously, under the same replacement of ai, 02 by Ai/7^, ^2/7!) respectively, we can take advantage of Eqs. 
(36) and (18) (with Ki = K2 = 0) to let Eq. (41) become, 

[A^A2{x'^y'^ + j2X^y + jixy^ + jij2xy) - A27i^(y + j2)Xy]D + Ai^^x^ + A2'yfy^ + 'yh^z = 0. (43) 

Set 71 = 72 = 1- The elimination of x and y between Eqs. (42-43) yields a degree 7 polynomial condition for D, 

V{D, z, Ai,A2, X) = Al (Ai + AXf [{z + Xf + Aiz] [{z + X + A2f + Ax{z + A2) ] - 
4^2(^1 + 4A) [ {z + Xf + Aiz] [AiAl + 5A1A2X + SAjX + V1A2'>? -h 4A^ -|- (-3A1A2 - 4^^ 
4AiA-|-4A2A-F4A2)^-4(^i -h2A2-FA)^2 -4^3]£)6 -I- ... - 64^(^i -|- ^2 + 20A -|- 16^) £) 256^ = , (44) 

which is too cumbersome to be listed here entirely. A factor \(z + A)^ + A\z\ forces its coefficients for both rP and 
to vanish when z = 770, see Eq. (40). Hence two roots D diverge at the expected one-body threshold. We also 
notice that for z = the two lowest degree coefficients of V vanish, hence a double root Z? = occurs. But, for 
the sake of simplicity in this Section, we shall not elaborate much on the exact nature of this two-body threshold 
singularity for this second model. Similarities with the behavior of the first model around 2; = are likely. In the 
following we rather study in some detail the singularity at z — r\Q. 

The degree 7 for D is familiar from the model of the previous Section. But the degree for z is now 5 rather than 
4. We verified that the limit A factorizes V{I), z, A\, A2, A) into a factor z and a polynomial with degree 4 for z. 

It is convenient to set special values for a numerical investigation, for instance ai = 71 = Ai = 72 = 1 and 
a2 = A2 = X = 2. The full polynomial then reads, 

r = 3{l + z) (A + z) [27 {3 + z) {6 + z) D + 12 (-108 -9z+ 14z^ + 2z^) ] - 
4 ( -2808 -I- 4770 -|- 59840^ + 3480^^ -|- 690^^ -|- Uz^ ) + 8(-584 -|- 1760z -|- 3323^^ + 1387^^ -|- 200^^ + 8z^)D^ - 
8 ( -140 -I- 11950 -I- 14710^ + 420z^ + 32z^ ) + 32(-4 -|- 107^ -|- 82z^ + 12z^)D^ - 16^(43 + 16z)D + QAz. (45) 

Here, the bound state lies at 770 = —1 with = The second solution, wq = —2i,r}o = —4, of Eqs. (39-40) violates 
the condition ^x > 0, hence pertains to an unphysical sheet. 

The 7 trajectories shown in Fig. 6 are those of the roots of V, Eq. (45), when runs from —7.5 to +4. This range 
suffices here to obtain a reasonable estimate of the root behavior when the energy runs from —00 to -l-oo. For the sake 
of graphical convenience, a renormalization D/{1 + \D\) forces large D's back to the trigonometric circle. Also a small 
imaginary part Qz = .2 is set to enforce the rule < 0. Black dots (or lines when nearing dots fuse) correspond 
to Siz < —4.1. Green and red ones correspond to —3.9 < Sftz < —1.1 and —0.9 < SRz, respectively. Finally, blue and 
yellow ones investigate neighborhoods, —4.1 < SRz < —3.9 and —1.1 < 5f 2: < —0.9, of expected singularities at z = — 4 
and z = —1, respectively. It turns out that Fig. 6 does not yield much information out of such "blue" and "yellow' 
segments, although it is clear that only two "loops" satisfy both rules 3D < 0, V3?z and lim^^.^^^ D = 0. Clearly, 
for a thorough investigation of all branchings and divergences, we should eliminate D between P and its derivative 
dV/dD, then study the neigborhoods of all the roots of the obtained resultant, 
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(42a) 
(42b) 

(42c) 



7^ = z (1 + zf{2 + z)(3 + z)(4 + z)2(6 + z)(32 + 7zf x 
(2426112 + 17293824Z + 540267842^ + 121209152z^ + 233641545z^ + 328920768z^ + 
307812074z^ + 191171112z^ + 79534245z* + 21923392z^ + 38269442^° + 380928z" + 163842^^ f . (46) 

This straightforward but lengthy task gives too cumbersome results to be published here, naturally. Still, it might 
be useful to compare Fig. 6 with a superposition of Figs. 2 and 11, keeping in mind that symmetry breaking double 
roots of the previous model will now be disentangled. Indeed, under the already mentioned two criteria, namely i) 
Z? ^ if |z| — > cx), and ii) 3D < 0, only the "tiny" loops selected from Figs. 2 and 11 survive. Letting 5z 0, we 
obtained graphical evidence that such two loops grow in such a way that their "blue" and "yellow" segments show 
the diverging roots predicted from the factor (1 + z)(4 + z) in front of D'' and . With the same renormalization 
D/(l + \D\), Fig. 7 confirms that two branches only are compatible with rules i) and ii), when we freeze "Slz — —1 
and let 3z > run. One of the "good" candidate roots diverges for z — —1, see the green segment in the lower left 
part of Fig. 7. The other "good" candidate, Di ~ .31 — .21i, see the small green segment at the beginning of the 
smallest trajectory in the lower right angle of Fig. 7, is a simple root as a function of z in this area, and deserves 
little comment. The diverging root, however, because of its quadratic branching, deserves a study of its reciprocal, 
d = D^^ . We set z = —1 + Z and expand V, Eq. (45), at lowest orders with respect to d and Z, 

V{d-^, Z - 1, = 1, As = 2, A = 2) = 270 {2d^ + 9Z) + 0{Z d). (47) 

The neglected term is of order Z"^/^, because, obviously, the leading order of the double root is d = ±3(— Z)^/^/\/2, 
real below and imaginary above threshold, respectively. A similar, straightforward argument for the vicinity of the 
additional, but unphysical threshold at z — —4 yields the leading order d = ±3(— 4 — z)^/^/(2\/2). 




FIG. 6. Complex D-plane. Trajectories of the 7 roots of Eq. (45) when —7.5 < SR^ < 4 and — .2. Black dots or 
lines correspond to SRz < —4.1. Blue, green, yellow and red ones correspond to —4.1 < SRz < —3.9, —3.9 < SRz < —1.1, 
— 1.1 < 5R2: < —0.9, and —0.9 < 5Rz, respectively. Only two trajectories always keep QD < and cancel D when \z\ 00. 
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FIG. 7. Complex D-plane. Trajectories of the 7 roots of Eq. (45) if = —1. Green, yellow and red mean < < .01, 
.01 < Qz < 2 and 2 < Qz < 8, respectively. Again, only two trajectories maintain < and cancel D when \z\ oo. 

Among all the singularities of this second model, we shall mainly discuss the physical threshold z = — 1. Eliminate 
z between Eq. (42a) and Eq. (42b), or, equivalently, subtract the equations, Eqs. (37), from each other, hence 

du}2 duji 

a2^2J2-r-r ^ aiUJiIi-f^, (48) 

a relation similar to Eq. (16). To prove the statement that the static HF energy 770 indeed defines a threshold solution 
of the TIME equations, Eqs. (37), it is enough to set 771 — > 770, 7/2^0 and z — > 779. This automatically induces 
z — aiUJi — a2^^2 ^ 0; naturally. Set 71 = 72 = 1, for a trivial scaling. Then Eq. (48) reads, 

M y^l + y) ^ X [Ai x{l + x)-\] 

l + 2y l + 2x ' ^ ' 

When Ai — 1 and ^2 = A = 2, we know that the limits of interest are a; — > 1, j/ ^ and z —1. These satisfy the 
condition, SRx > 0, hence only diy must be investigated. Define X = x — 1 and Z = z + 1. Then Eq. (49) boils down to 
2y^ = X, at leading orders in y and X. Accordingly, Eq. (42b), for instance, boils down to, 2y^ + Z = 0. For z < — 1 
(below threshold) the solution, y — > y/ {—z — 1) /2, is acceptable, with, simultaneously, x —z. For z > ~1 (above 
threshold), however, we find that a small, but positive Sjz is necessary to allow the condition JRy > 0. This occurs 
because for Z > the leading order, y^ —Z/2, actually generates 3?/ only. An expansion up to higher orders, is 
thus necessary for the knowledge of SRy. A straightforward, but slightly lengthy calculation yields, 

y = zZ'- '-Z'^ - h'^ + + O , (50) 



where Z', a positive number, is defined as Z' = ^ Z/2 = yj [z + l)/2. The "formal conjugate" of this expansion, 

y = ~^Z' + '-Z'^ - + O , (51) 

also holds, naturally. (Equivalently, it means the opposite choice of Z' , namely Z' = Z/2.) Both expansions 
induce a negative 3??/ as long as Z is real and positive. The sign of this ^y can be easily reversed, however, as soon 
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as, above that threshold z = —1, an imaginary part is implemented. Another slightly cumbersome calculation 
defines, upon taking advantage of either Eq. (50) or Eq. (51), the condition for the border at which one of such 
roots acquires a positive real part. This is illustrated by Fig. 8. Near to that threshold z — —1, the leading orders 
of the border condition give, (3z)^ = |(3?2; + 1)^. Other numerical values for the parameters 7i, etc. modify the 
numerical analysis, naturally, but leave intact the conclusion, namely that \^z\ must have at least a non vanishing 
value above the threshold if one needs one of these two roots to be compatible with the condition, 3?y > 0. 



.14 



. 12 



.01 



0.06 



. 02 








FIG. 8. Complex Z = z + 1-plane. Below the plotted line, both roots described by Eqs. (50-51) show "^y < 0. Above that 
line, one of them shows 5Ry > 0. The line contains both thresholds Z = and Z = 1. 
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FIG. 9. Complex y-plane. Trajectories of y when 5Rz = — .6 and increases from 0. Blue lines are trajectories for which 
either Ka; or 5Rj/ or both are negative. Only one branch, that long one in thelower right quadrant, survives the double condition, 
Ka: > 0, 5Ry > 0. Tiny blue segment, < 9z < .03. Green segment, .04 < < .2. Orange one, .3 < Qz < 1.8. Red one, 
2 < 4. 
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We show on Figure 9 the trajectories of y for < 3z < oo when = —.6 is frozen at an intermediate value 
between the thresholds z = —1 and z = 0. Only one branch is of interest, because all the other branches either stay 
in the 5iy < sector or the x partner root shows 3?a; < 0. The tiny blue segment at the beginning of this branch 
corresponds to < .03, imaginary parts too small for letting y acquire a positive real part, see Figure 8. 



IV. A THEOREM 



We return to the case where N is any finite particle number. The two-body interaction V = J2i>j ^ij contained 
in the physical Hamiltonian H is assumed to be made of short ranged potentials Vij . Then the TIMF mean fields Ui 
are also short ranged. For details of a further antisymmetrization with identical fermions, where the mean potential 
will be the same U for all particles, we refer to [8]; the short range of U remains, whether one considers its direct or 
exchange part. At present we still retain the case of distinct particles. Eqs. (3) read again, 

iVi-ti-Ui)\^i) = \Xi), {<fi\iVi-ti-Ui) = {x'i\, (52) 

with 

i^'m) {^',\{ti+Ui)\^i) 

Notice that we now process a generalized argument, since we can also study non diagonal elements {x'\{z — H)~^\x) 
where x a-nd x' a-re products made of orbitals Xi ^nd x'i, respectively. The Euclidian restriction is not implemented 
any more. The trial functions </> and </)' are the products made of orbitals (pi and ip'^, respectively. All such quantities 
and wave functions depend on z, but we stress here that, because of the short range of Ui, the spectrum of hi = ti + Ui 
has a fixed continuum, extending from to +oo on the real axis of the rji complex plane. In general Ui is complex 
and the poles of {r]i — hi)~^ need not be real; as a matter of fact they move as functions of z (and of the choices of x 
and x')- But the continuum cut for the spectrum of hi remains always the same. It is therefore legitimate to ask the 
question "what happens if one of the rjiS vanishes, hitting the threshold of the continuum of hiT\ Incidentally, it will 
be noticed that there are many trajectories (sheets) of such ry^'s as functions of z. The multiplicity comes not only 
from the existence of N "momenta", Wj oc ±v^-, with their ± ambiguity [12], but it is also due to the nonlinearity of 
the mean field theory. For instance, in our second model, we found seven sheets, see the seven roots for each quantity 
D{z), x{z), y{z) driven by z. 

As a preliminary remark, we use Eqs. (53) to notice that the mismatch between any propagation energy r/j and the 
corresponding self energy {fi\{ti + Ui)\(pi) / {f'i\ipi) does not depend on i. Furthermore we can take advantage of Eqs. 
(52) to relate the self and propagation energies as, 

= + K:i = {x^lihi - ,ir\xi) = -(^^Ix.) = -bdl^i) . (54) 

In other terms, the mismatch is measured by the ratio {x!i\'fi) I {'f'iWi) = {'f'i\Xi) / ^ a- function of r]i. When 
calculated at self consistent r]i{zys, such ratios do not depend on i any more. 

Assume that the special vanishing rjs reads rjg — ie^ , where e is a real, positive infinitesimal. This means that 
we select in the z complex plane a trajectory which in turn induces an 7/^ trajectory leading to retarded, outgoing 
boundary conditions for that special p>s and its partner t^^. For the sake of simplicity, set the inverse mass coefficient 
as to unity, or, equivalently, renormalize rjs and Us accordingly. In physical three dimensions, the partial wave 
components ipsi are described by differential equations of the form. 



1) 



+ Useir) - is 



/•oo 

^se{r) = Xse{r) - V / dr' Usee' {r, r') psi' (r'), , (55) 

e, Jo 



where Usee' is a short notation accounting for, if necessary, partial wave coupling and/or non local parts of Us- The 
source term Xs is expanded in partial waves as well, naturally. It is then convenient to denote the right hand sides of 
Eqs. (55) as source terms ^seir)- These are short ranged, obviously again. Similar equations hold for <p^. 

For each £ let (7se{'>') be the regular solution, usually normalized as (r'^fXO) = 1, of the homogeneous, left hand side 
of Eqs. (55). The short range of Ug, and similar short ranges assumed for x and x', make it that, when r oo, then 
fseir) becomes ~ exp[(i — l)er/\/2] dr' ase{r')^se{'>''), with a similar asymptotic formula for tp'^^. Let C, a real 
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and strictly positive number, be any convenient lower bound for the absolute values of these integrals J and J 
in a neighborhood of e ^ 0. This C exists, since such integrals are usually finite and non vanishing when e = 0. It is 
clear that, as e — > 0, there are no more any exponential decays or any asymptotic oscillations in the product 
Then, at this limit for e, the integral (</?g|<^s) diverges, while obviously an integral such as (Xsl'^s) remains finite. The 
"mismatch" cancels out. 

This indicates that, for any i ^ s, the ratios ICidrji/dlCi vanish simultaneously at their respective energies 77^. 
Besides threshold limits for each r^i, there is an easy interpretation for such a situation, namely, each among such 
N — 1 propagation energies converges towards a bound state energy of its hi. Indeed, let dei be an infinitesimal 
difference between r/i and an isolated eigenvalue of hi. Then it is trivial, in an energy representation with biorthogonal 
eigenstates of /i^, to sec that ICi diverges at order {dei)~^ ^ while dJCi/drji diverges at order (dci)"^. 

The situation is thus representative of a Hartree(-Fock) solution for the TV — 1 particle system. This is confirmed by 
the observation that, since {ip'g\(fis) diverges, the potential f/j/g induced by particle s upon any particle s vanishes. 
Indeed, the short range of v in the formula, 

u,sin)= ^'''-''-^'^^7::^f''^^^^'-'\ (56) 

makes the numerator converge Vrj, while the denominator diverges. Any matrix element {(p^(p'g\v\(pi(ps) will vanish 
too, for the same reason. Furthermore, the full matrix element {(j)'\H\(f)) / {(j)'\(p) can always be split as. 



m {Ct>'-s\<l>-s) {v'sWs) 

where the subscript —s refers to the subsystem where particle s is removed. At the limit under study, both rjs and 
ICs drjs/dlCs vanish. Honco, according to Eq.(54), the self energy for particle s vanishes and the full matrix clement 
{(j)'\H\(f)) / {(f)'\(f)) reduces to the subsystem value, {(j)'_g\H-s\(j)-s) / {4''-s\'t'-s)- Furthermore, setting i = s in Eq.(53), we 
find that z {(j)'_g\H-s\(l)-s) / {4>'-s\4'-s) ■ The threshold for the continuum of particle s in the 2-plane corresponds to 
the Hartree(-Fock) binding energy of the subsystem. 

Conversely, if z converges towards a Hartree(-Fock) bound state energy of an iV — 1 particle system, it is easy to 
verify that at least one solution of the TIMF equations for the N particle system consists in a threshold wave for the 
additional particle, as a spectator of the static solution for the subsystem. 

Notice that several special particles, not just one, can be forced into their continuum thresholds simultaneously. 
For instance, if particle s and s' are such that r/s = rjs' = 0, then all potentials Ui/g and Ui/s', including Ug/s' and 
Ug'/s, vanish, and z = {(f>'_g_g,\H-s-s'\(l>-s-s'} / {(l>'-s-s'\'l^-s-s'} , a subsystem energy for A'' — 2 particles. 

It can be also noticed that such singularities do not depend upon the source terms x and Indeed, the locations 
of such thresholds derive from homogeneous equations, where only H appears. 

The present theorem can be phrased in a way which generalizes the theorem of [10]: not only the mean field binding 
energies of a system of N particles define singularities of the TIMF propagator, but the mean field binding energies 
of its subsystems define thresholds of cuts where the additional particles become unbound. 



V. DISCUSSION AND CONCLUSION 



There are two parts in this work, namely on the one hand a couple of very special, analytical models, see Sections 
II and III, and on the other hand a theorem of a more general validity. 

The systems described by our models are physically trivial, since they make non interacting particles. But their 
mathematical interest is different. As stated at the beginning of this work, it is important, for large particle numbers, 
to validate the replacement of convolutions by straight products, and our models allow a detail study of all singularities 
and nonlinearities introduced by the mean field approximation. We investigated three representations, namely what 
happens in, i) the z-plane (propagation energy), see for instance Figure 4, ii) the D-plane (TIMF amplitude), see for 
instance Figure 1, iii) pseudo momentum planes, such as, for instance the case of x = —i\/r]i/ai, see Figure 5. Since 
our models automatically implement an analytic continuation from physical to unphysical sheets, there is no cut to 
consider in the pseudo momentum complex planes. It is obvious, however, that for both pseudo momenta x and y 
the imaginary axis represents both rims of the cut which would be necessary in their respective Tj-plane. Accordingly, 
see for instance Figure 8, values of z for which the real part of a pseudo momentum vanishes, or identically for which 
a propagation energy r]{z) becomes real and positive, make cuts in the z representation. The zoology of the TIMF 
solutions turns out to be surprisingly rich. The main two conclusions provided by the models can be listed as follows. 



16 



i) except when the many-body propagation energy z has too small an imaginary part, the TIMF equations always 
generate at least one branch of solutions where each single particle undergoes a retarded propagation and the TIMF 
amplitude D shows all suitable properties needed for a reasonable approximation of a Green's function matrix element, 

ii) the threshold of a single particle continuum induces the threshold of a cut singularity in the z representation; if 
one calls "projectile" that special particle becoming unbound, and "target" the system made by the other particle, 
the corresponding threshold value for the full propagation energy z is the binding energy of the "target" . 

The theorem derived in Section IV, valid for any particle number N > 2, extends this numerical and analytical 
evidence. Hence the mean field theory of collisions mimicks the connection between singularities of the inhomogeneous 
problem {z — H)\^) = \x) and the solutions of the homogeneous Schrodinger equation {E — H)\^) = 0. At this stage of 
our work, the similarity is restricted, however: we considered only partitions where a "target" is surrounded by one or 
several unbound particles, and we have not proven thresholds defined by mean field energies of partitions N1+N2 = N, 
Ni > 2, N2 > 2, into two clusters, each of them carrying its full internal energy. Nor have we considered even finer 
partitions Ni + N2 + N3 = N, with A^i > 2, N2 > 2, N3 > 2, and so on. Last but not least, the present work lacks 
a clear description of the shapes of the cuts beyond their thresholds. The preliminary resTilt obtained at the stage of 
Figure 8, with a "border equation" like {Qz)"^ = |(3?z + 1)^, is an omen of subtle arguments yet to be phrased. 

Despite such questions still open, the TIMF approximation now appears like a theory of collisions endowed with 

properties, such as poles and thresholds, with soimd interpretations in terms of Hartree-(Fock) energies of subsystems. 
The special role played by single particle energy propagators {rj — h)^^ in the definition of such properties is a logical 
consequence of the factorization of trial wave functions, an essential ingredient of practical approximations. With 
the present and foregoing studies, TIMF appears as a reliable and practicable alternative to resonating group (RGM) 
or generator coordinate (GCM) studies for application in nuclear astrophysics where there is still a demand for 
microscopic rather than phenomenological calculations of processes relevant to element synthesis. 
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For the symmetry breaking sector of the first model, we again set fli = 7i = 1, scaling x, y, z and ensuring the 
factorization of the resultant, Eq. (22a). It is easy to analyze the singularities of the direct solution of Eq. (26a), 
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(58) 
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in terms of one cut from —1 to in the complex z-plane, or, alternately, two cuts from — oo to —1 and from to +oo, 
and observe that the square root singularity at z = seems to represent a very traditional threshold singularity. Less 
physical, the role of = —1 is to reflect the discriminant A2 = ^zj (1 + z) of Eq. (26a). 

In the forthcoming Figures, we keep = 1, as a natural scale for energies and inverse amplitudes. Fig. 10 shows 
the graph of z^/c when IJ^fc is real and takes on all values from —00 to +00. The symmetry axis at = 2 is obvious 
from Eq. (27). Since the physical amplitude is negative when z is negative, the right lower branch of the graph is 
clearly unphysical, while the left lower branch is a reasonable candidate for approximations. (Notice, however, that 
no real estimate of the amplitude is offered for — 1 < 2; < 0.) In turn, the right upper branch is also ruled out, as T) 
must vanish when z — > +cx). This leaves the left upper branch as a tolerable candidate for physical approximates of 
the real (principal) part of "D when z is positive. 




FIG. 10. Scaled energy z (unit &) as a function of the "symmetry breaking" amplitude D (unit 1/^). 

Rather than considering inverse functions z{D), we then show in Figs. 11-12 the trajectories, in a complex plane 
"_D", of the solutions of Eq. (26a) when takes on all values from —00 to +00 and 3z is frozen at some fixed 
value F. The physical situation corresponds to F = 0+, naturally, but Figs. 11-12 use larger values of F for graphical 
convenience. For Fig. 11, we use F = .075 and F — .4, which generate for Df,fc two "outer loops" and two "inner 
loops", respectively. The role of D = 2 as a symmetry center is obvious. The shrinking of the loops when F increases 
comes from the fact that, as \z\ 00, the dominant part of the symmetry breaking equation is D{D — 4) = 0. 
Conversely, the evolution of such loops into "angles" when F — > is transparent on Fig. 12, obtained with F = .02. 
Only those solutions which lie in the lower half plane can be retained as physical candidates, according to the condition 
"if 3z > 0, then 31? < 0." Hence the general physical behavior of Dbk is as follows: 

- when z is real and increases from —00 to —1, then D decreases from 0^ to —00, 

- when z is real and increases from —1 to 0, then D varies from 2 — ioo to 2, 

- when Qz = 0+ and 3iz increases from to -|-oo, then D decreases from 2 + iO~ to iO~ . This infinitesimal imaginary 
part hints that D^k can at best approximate the principal part of V. This was already deduced from Fig. 10. 
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FIG. 11. Complex D-plane. Trajectories of symmetry breaking D when 9z = .075 (outer loops) and "Qz = .4 (inner ones). 
For S$2 = .075, blue dots growing for growing between — oo and 0^ and red ones growing for growing from 0"'" to +oo. 
For Ssz = A, blue and red replaced by green and yellow, respectively. 
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FIG. 12. Complex D-plane. Lower loop trajectory of Dhk when Sjz = .02 Blue dots growing for 5R2: growing between —00 
and 0~. Red ones growing for SRz growing from 0^ to +00. 

The "breaking" factor of the factorizing resultant between Eqs. (20) reads, 



19 



x^-2x^z-x'^z-2xz{l + z)-z{l + z)=0, 



(59) 



keeping in mind that its roots must be paired as {x,y). Fig. 13 displays a contour plot of the product of the 
corresponding 4 real parts of the roots as functions of z. The corresponding cut in the 2;-plane is the border between 
the light grey and the darker grey areas. It is now made of two branches. The right hand branch, while not located 
on the real axis of the ^-plane, again contains the two-body threshold z = 0. 
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FIG. 13. Complex z-plane. Cuts caused by the condition SRx = for symmetry breaking roots. The cuts are the borders 
between light grey and darker grey areas. 

Then Fig. 14 shows the trajectories of the 4 roots when we freeze 5ftz = 0.1 and let 3z nm from —1 to +1, allowing 
z to cross twice the right hand side cut shown by Fig. 13. The sizes of the dots are coded like those of Fig. 5: minimal 
for ^z = —1, growing until '<sz = 0, minimal again for small positive values of growing again until 'Qz = 1. The 
"vertical" branch on the right hand side of Fig. 14 has the unsatisfactory property, that its "t/-partncr" according 
to Eq. (22b), is the loop like, tiny branch on the left hand side of Fig. 14. Hence Uxdiy < 0. In turn, the two 
"horizontal" branches on Fig. 14 are "x-y" partners and are partly located inside the right hand side of the complex 
X, y-plane. But actually the double condition, ^x > 0, 5fj/ > 0, is never satisfied. It must be concluded that the 
symmetry breaking sector is unphysical. 
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FIG. 14. Complex x-plane. Trajectories of the symmetry breaking a;'s when = 0.1, while Qz runs through one of the cuts 
shown by Fig. 13. Blue dots growing when grows from —1 to 0. Red dots growing when Qz grows from to 1. 

A trivial manipulation of Eqs. (20) shows that the pairing of roots, for this sector and such special parameters, 
follows the rule. 



xy 



X — y, hence y — 



x + 1' 



(60) 



which is its own inverse tranform, naturally. The rule is equivalent to Eq. (22b), but simpler. Then if one defines 
s = X + y, it is easy to reduce Eq. (59) into. 



(z-s)2 = z(z + l), or z-s = ±[z(z-Hl)]i/^ 

while shortening Eq. (60) into, 

xy^z-s^ ±[z(z + l)]^/2. 
This means that Eq. (59) factorizes into two distinct equations, 

X^^x\z-{Z + z2)l/2l + + 22)1/2 ^ ^ 
Z + (z + Z^)'/^] ~ [Z + z')'^^ = . 



(61) 
(62) 

(63a) 
(63b) 



It is easy to verify that each of these is invariant under the transform, Eq. (60), hence each yields a pair {x,y). A 
detailed analysis of all cases for such equations is trivial, but too lengthy to be published. Rather, it is enough and 
easy, actually, to set 3z = 0, and plot, for instance for Eq. (63a), its two numerical roots as functions of z. It turns 
out that at least one of the roots has always a negative real part. The same phenomenon occurs for Eq. (63b). All 
told, the symmetry breaking sector does not respect the constraints requested simultaneously for ^x and 5fty. 
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